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ABSTRACT:  Interest  in  materials  that  undergo  singlet  fission  Singlet  Fission: 

(SF)  has  been  catalyzed  by  the  potential  to  exceed  the  A  Direct  Mechanism 

Shockley— Queisser  limit  of  solar  power  conversion  efficiency. 

In  conventional  materials,  the  mechanism  of  SF  is  an 
intermolecular  process  (xSF),  which  is  mediated  by  charge 
transfer  (CT)  states  and  depends  sensitively  on  crystal  packing 
or  molecular  collisions.  In  contrast,  recently  reported 
covalently  coupled  pentacenes  yield  ~2  triplets  per  photon 
absorbed  in  individual  molecules:  the  hallmark  of  intra¬ 
molecular  singlet  fission  (iSF).  However,  the  mechanism  of 
iSF  is  unclear.  Here,  using  multireference  electronic  structure 
calculations  and  transient  absorption  spectroscopy,  we 
establish  that  iSF  can  occur  via  a  direct  coupling  mechanism 

that  is  independent  of  CT  states.  We  show  that  a  near-degeneracy  in  electronic  state  energies  induced  by  vibronic  coupling  to 
intramolecular  modes  of  the  covalent  dimer  allows  for  strong  mixing  between  the  correlated  triplet  pair  state  and  the  local 
excitonic  state,  despite  weak  direct  coupling. 


Singlet  fission,  the  process  by  which  one  singlet  exciton  splits 
into  two  triplet  excitons,  is  proving  an  important  property  for 
materials  used  in  third-generation  solar  cells  and  photodetectors, 
among  other  optoelectronic  devices.1  Unfortunately,  techno¬ 
logical  applications  are  limited  by  the  small  number  of  organic 
chromophores  that  undergo  efficient  singlet  fission,4’5  and  are 
largely  restricted  to  molecular  crystals  of  oligoacenes  and  related 
materials  that  exhibit  intermolecular  singlet  fission  (xSF).6-10 
Insights  into  the  mechanism  of  xSF  have  been  obtained  by  tuning 
spatial  interchromophore  interactions;6’11  14  however,  the 
sensitivity  of  the  mechanism  to  crystal  packing  orientations 
poses  a  significant  challenge.  As  such,  devices  based  on  xSF  are 
limited  by  the  lack  of  high  throughput  processing  strategies  to 
create  highly  ordered  molecular  structures.  In  contrast,  recently 
reported  efficient  intramolecular  singlet  fission  (iSF)  materi¬ 
als15  "2  offer  significant  advantages  in  terms  of  their  tunable 
molecular  and  electronic  structure,  solution  processability,  and 
the  ability  to  form  tailored  interfaces.  ~  However,  a  limited 
understanding  of  the  fundamental  process  of  triplet  pair  (2XT:) 
formation  hinders  the  design  of  versatile  iSF  materials.  It  is 
therefore  imperative  to  uncover  the  detailed  mechanism  of  iSF 
and  to  establish  its  relationship  to  chemical  structure.  Current 
literature  on  the  mechanism  of  xSF  is  in  general  agreement  that 
charge-transfer  (CT)  states  play  a  significant  role  in  mediating 


the  coupling  between  local  excitonic  (LE)  singlet  states  formed 
upon  photon  absorption  and  the  correlated  triplet  pair](TT)  or 
multiexcitonic  (ME)  state."  In  addition,  the  CT-mediated 
mechanism  has  been  implicated  in  iSF  polymers,  using  strong- 
donor/strong-acceptor  interactions  along  the  conjugated  back¬ 
bone  and  leading  to  high  yields  of  triplet  pairs  upon 
photoexcitation  in  dilute  solution.16,  However,  iSF  materials 
that  operate  by  this  mechanism  can  suffer  from  relatively  short 
triplet  lifetimes  and  parasitic  loss  pathways  involving  charge- 
separated  and  dark  singlet  states,16  decreasing  the  net  yield  of 
triplet  pairs. 

Recently,  covalently  coupled  pentacene  dimers  were  synthe¬ 
sized  and  demonstrated  to  yield  triplet  pairs  nearly  quantitatively 
(~200%  triplet  yield)  by  ultrafast,  sub-picosecond  iSF.  '  ’ 22 
From  a  practical  standpoint,  bipentacenes  are  the  first  family  of 
iSF  chromophores  with  excited  state  dynamics  on  time  scales 
comparable  to  xSF  in  crystals.  Although  it  is  known  that 
pentacenes  meet  the  energetic  criteria  for  singlet  fission,  where 
the  energy  of  the  singlet  state  is  at  least  twice  the  energy  of  the 
triplet  state,  EfSj]  >  2E[T,],  the  mechanism  of  iSF  in  pentacene 
dimers  is  unclear.  In  particular,  the  role  of  CT  states  in  mediating 
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iSF  is  called  into  question:  low-lying  CT  states  are  not  expected 
based  on  the  molecular  structure,  and  experimental  measure¬ 
ments  of  iSF  rates  show  relative  insensitivity  to  changes  in 
solvent  properties.  This  raises  the  possibility  of  an  alternate  direct 
coupling  mechanism,  where  the  correlated  triplet  state  is  directly 
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Figure  1.  (A)  Direct  (in  green)  and  charge  transfer  mediated  (in  red) 
mechanisms  of  singlet  fission.  (B)  Molecules  studied  in  this  manuscript, 
where  the  abbreviations  represent  bipentacene  and  its  dihedral  angle  9 
(BP-0)  where  9  =  37°,  42°,  and  57°  for  the  various  twist  angles.  (C) 
Schematic  excited  state  rates  that  depend  on  dihedral  angle,  8,  after 
excitation  (hi/)  to  the  singlet  state  (S, ):  rate  of  intramolecular  singlet 
fission  (rlSF)  and  the  rate  of  triplet  pair  (2xT,)  recombination  (rrec). 


rates  of  iSF  and  triplet— triplet  recombination  vary  dramatically 
as  a  function  of  the  coupling  interactions  in  the  bipentacene,  it  is 
clear  that  structure— property  relationships  can  be  exploited  to 
control  excited  state  dynamics  and  to  test  different  mechanisms 
for  the  SF  process.’3 

Here,  we  start  by  investigating  the  direct  coupling  mechanism, 
where  a  vibronically  induced  degeneracy  between  the  local 
excitonic  (LE)  and  multiexcitonic  (ME)  states  leads  to  ultrafast 
iSF  despite  a  small  coupling  matrix  element  between  the  LE  and 
ME  states.  We  identify  a  specific,  high-frequency,  intramolecular 
vibrational  mode  that  assists  direct  population  transfer  from  the 
LE  state  to  the  ME  state  through  an  avoided  crossing  (see 
Supporting  Information  movie).  This  demonstrates  that  ultrafast 


iSF  can  occur  via  the  direct  coupling  mechanism  in  covalently 
linked  pentacene  dimers  and  does  not  require  CT  mediation.  We 
note  that  an  avoided  crossing  previously  identified  along  an 
intermonomer  mode  in  cofacial  pentacene  dimers  was  shown  to 
be  an  artifact  of  the  calculation  method,32’36  and  although 
intramonomer  ring-breathing  modes  have  been  identified  as 
important  in  xSF,3  1,40  their  role  has  not  been  clearly  established. 
We  also  investigate  the  possibility  of  a  CT-mediated  super¬ 
exchange  mechanism  for  iSF  in  bipentacenes.  We  show  that 
although  iSF  rates  sensitively  depend  on  the  dihedral  angle,  CT 
state  energies  are  largely  independent  of  this  factor.  Experimental 
measurements  establish  the  independence  of  iSF  rates  on  solvent 
properties,  specifically  orientational  polarizability.  Furthermore, 
using  both  theory  and  experiment,  we  quantify  the  role  of 
structural  factors  (Figure  IB)  by  studying  the  change  in  rates  of 
iSF  with  the  dihedral  angle,  6,  between  monomers  (Figure  1C). 
While  the  direct  coupling  mechanism  correctly  predicts  the 
observed  change  in  iSF  rate  with  dihedral  angle,  the  indirect 
couplings  that  drive  CT-mediated  singlet  fission  do  not  show 
concomitant  change.  Taken  together,  these  insights  into  the 
mechanism  of  iSF  from  experiment  and  theory  suggest  a  new 
class  of  materials  where  the  direct  coupling  mechanism  is 
important. 

■  THE  EXCITED  ADIABATIC  AND  DIABATIC  STATES 
OF  BP-37 

We  begin  our  theoretical  investigation  with  the  parent 
bipentacene,  BP-37,  where  the  notation  indicates  the  calculated 
dihedral  angle  of  the  minimum  energy  structure,  9  =  37°.  This 
compound  is  a  simplified  version  of  the  experimentally  measured 
BP-2H,  which  undergoes  iSF  with  a  0.76  ps  time  constant  in 
~200%  triplet  yield,  followed  by  triplet  pair  recombination  with  a 
time  constant  of  450  ps. 1  In  the  truncated  version,  the 
solubilizing  tri-isopropylsilyl  groups  are  replaced  with  hydrogens. 
We  establish  that  the  solubilizing  groups  present  in  the 
experimental  study  do  not  affect  the  electronic  structure  of  the 
monomer  significantly  (Supporting  Information),  in  agreement 
with  previous  work.’1  We  optimize  the  geometry  of  the 
substituted  pentacene  monomer  and  establish  that  its  geometry 
changes  little  upon  dimerization.  This  allows  us  to  define  the 
dimer  reference  geometry  as  two  substituted  pentacene 
monomers  linked  through  a  covalent  bond  of  length  1.48  A  at 
the  2-  and  2'-positions  (Figure  IB).  The  excited  state  energies  of 
BP-37  are  obtained  from  complete  active  space  self-consistent 
field  theory  (CASSCF)  calculations  with  extended  multi¬ 
configuration  quasi-degenerate  perturbation  theory 
(XMCQDPT).3-’4  The  adiabatic  eigenstates  obtained  from 
these  calculations  are  labeled  S„  where  i  =  0  indicates  the  ground 
state  and  higher  values  of  i  represent  the  excited  states. 
Throughout  this  manuscript,  we  also  employ  quasi-diabatic 
states  constructed  by  rotating  the  adiabatic  states  into  a  local  state 
basis.  The  reference  diabatic  states  are  labeled  following  standard 
notation  where  the  first  index  of  the  label  refers  to  monomer  A 
and  the  second  index  refers  to  monomer  B  in  the  dimer.  The  low- 
lying  diabatic  states  are  shown  in  Figure  2  and  correspond  to  the 
ground  state  gg,  the  ME  state,  tt,  the  two  singlet  LE  states,  eg  and 
ge,  the  two  CT  states,  ac  and  ca,  and  the  doubly  excited  states,  dg 
and  gd.  We  will  use  these  adiabatic  and  diabatic  states  to  aid  in  the 
analysis  of  both  the  direct  coupling  and  CT-mediated  super¬ 
exchange  pathways  for  iSF  in  bipentacenes. 
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Figure  2.  Eight  lowest  diabatic  configurations  (excited  states  in  gray 
boxes)  represented  using  the  highest  occupied  molecular  orbital 
(HOMO)  and  lowest  unoccupied  molecular  orbital  (LUMO)  on  each 
monomer  (A  and  B),  highlighting  the  various  states:  charge  transfer 
(CT),  multiexcitonic  (ME),  doubly  excited  (DE),  and  local  excitonic 
(LE). 

■  THE  DIRECT  COUPLING  MECHANISM 

We  find  that  an  eight-orbital/eight-electron  (8o8e)  active  space 
is  necessary  to  yield  accurate  excited  adiabatic  state  energies  for 
the  dimer,17’4'1’44  and  our  results  are  shown  in  Table  1.  In  order  to 


Table  1.  Lowest  Vertical  Excitation  Energies  of  the  Covalent 

Dimer  at  the  Reference  Geometry 

state  8o8e  vertical  energy  (eV) 

S, 

1.96" 

s2 

1.96 

S3 

2.01 

aExperimental  bright  state  at  1.88  eV. 

quantify  the  direct  coupling  mechanism  we  perform  a  rotation  of 
the  adiabatic  states  into  a  localized,  diabatic  state  basis.29’30’45-4 
Unfortunately,  using  the  4-fold  diabatization  scheme  '  '  ' 
(Supporting  Information)  in  the  8o8e  active  space  yields 
inaccurate  quasi-diabatic  state  configurations  that  account  for 
<85%  of  the  higher  lying  adiabatic  states.  Therefore,  we  construct 
an  accurate  diabatic  state  representation  of  the  Hamiltonian  at 
the  8o8e  level  of  theory  by  including  only  the  three  lowest  lying 
adiabatic  excited  states  shown  in  Table  2.  We  analyze  the 


Table  2.  Diabatic  4-State  Hamiltonian  Matrix  Calculated 
Using  an  8o8e  Active  Space  (in  meV)  at  the  Reference 
Geometry  of  BP-37 


gg 

ge 

eg 

tt 

gg 

0.0 

-18.7 

-18.7 

-27.0 

ge 

-18.7 

1961.2 

0.8 

-2.2 

eg 

-18.7 

0.8 

1961.2 

-2.2 

tt 

-27.0 

-2.2 

-2.2 

2009.7 

contribution  to  the  excited  adiabatic  states  from  individual  state 
configurations  and  find  no  significant  CT  character  in  states 
within  1  eV  of  the  first  excited  state.  The  adiabatic  state  with 
primarily  ME  character  (Sj)  is  slightly  higher  in  energy  than  the 
adiabatic  states  that  have  primarily  LE  character  (S2  and  S3).  We 
find  that  there  is  weak  direct  coupling  between  the  LE  and  ME 
states, 


l(fflHlge)l  =  l(fflHleg)l  =  2.2  meV 

at  the  reference  dimer  geometry  and  that  the  energy  difference 
between  the  LE  and  ME  states  (~50  meV)  is  small,  as  compared 
to  previously  reported  values  for  other  acenes.12,30-32  The  near 
coincidence  of  energies  of  the  LE  and  ME  states  and  the  large 
energy  gap  between  these  states  and  the  states  with  significant 
CT  character  suggests  a  difference  in  the  mechanism  of  iSF  in  the 
present  molecules  from  an  earlier  study  of  xSF  in  a  staggered 
cofacial  pentacene  dimer.5" 

The  small  energy  separation  suggests  that  vibrational 
relaxation  might  induce  a  change  in  the  character  of  the  lowest 
singlet  state,  necessary  for  efficient  iSF.  The  next  step  is  to 
identify  key  vibrational  modes  that  couple  strongly  to  the 
electronic  states  and  to  evaluate  their  role  in  promoting  mixing  of 
the  LE  and  ME  states.  Since  BP-37  is  a  weakly  coupled  dimer,  we 
focus  on  the  normal  modes  of  the  substituted  pentacene 
monomer.  We  project  the  geometric  distortion  accompanying 
relaxation,  from  the  ground  state  structure  (vertical  excitation)  to 
the  optimized  S3  and  T,  monomer  structures,  onto  the 
corresponding  set  of  normal  modes.  We  find  a  dominant 
monomer  ring-breathing  mode  at  1435  cm-1  based  on  the 
magnitude  of  energy  changes  along  the  relaxation  pathway  and  a 
significant  Huang— Rhys  factor  (Supporting  Information).  Not 
surprisingly,  this  mode  is  consistent  with  the  vibronic 
progression  spectrum  of  an  unsubstituted  pentacene  mono¬ 
mer70’51  and  in  agreement  with  previous  work  that  suggests  the 
importance  of  a  ring-breathing  mode,39’40  although  an  avoided 
crossing  has  not  been  identified. 

We  use  this  ring-breathing  mode  to  generate  a  two- 
dimensional  potential  energy  surface  (PES)  for  the  excited 
states  along  the  dimensionless  coordinate,  Q,  corresponding  to 
the  two  monomer  normal  modes  at  1435  cm-1.  We  use  an 
inexpensive  and  efficient  occupation  restricted  multiple  active 
space  spin  flip  (ORMAS-SF)  approach  to  rapidly  scan  the  PES 
and  identify  geometries  of  interest.5"’  A  limited  number  of 
XMCQDPT  calculations  are  then  performed  with  an  8o8e  active 
space  at  these  geometries  to  obtain  accurate  excited  state 
energies.  Performing  a  diabatization  as  described  earlier,  we 
obtain  three  excited  diabatic  state  potential  energy  surfaces  along 
the  diagonal  of  vector  Q,  where  Q_  =  (0,0)  corresponds  to  the 
vertical  excitation  geometry.  Figure  3A  shows  that  as  we  move 
along  the  diagonal,  the  initially  higher  lying  ME  state  rapidly 
becomes  degenerate  with  the  LE  states,  and  then  significantly 
lower  in  energy.  At  the  point  of  degeneracy  (an  avoided  crossing) 
there  is  mixing  between  the  ME  and  LE  diabatic  states  that  leads 
to  a  change  in  character  of  the  lowest  adiabatic  excited  state  (Sj): 
Figure  3B  shows  the  composition  of  the  three  lowest  excited 
adiabatic  states  as  a  function  of  the  normal  mode  coordinates 
demonstrating  the  change  in  character  of  the  lowest  adiabatic 
state  as  a  function  of  Q. 

It  is  clear  that  despite  the  overall  weak  direct-coupling 
strength,  strong  mixing  occurs  as  the  adiabatic  states  become 
sufficiently  close  in  energy.  We  note  that  the  involvement  of  a 
high-frequency  mode  is  in  keeping  with  the  observed  ultrafast 
time  scale  for  iSF.  In  Figure  3C,  we  provide  physical  insight  into 
the  direct  coupling  mechanism  for  iSF,  mediated  by  a  vibrational 
mode.  Starting  from  a  vertical  excitation  at  Q_=  (0,0),  where  the 
bond  character  shown  is  deduced  from  the  reported  X-ray 
structure15  (Supporting  Information),  the  vibrational  mode  of  a 
substituted  pentacene  monomer  points  to  a  structure  where  the 
vectors  guide  the  shortening  of  the  carbon— carbon  bonds 
highlighted  in  red  (double-bond  character  in  Figure  3C), 
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Figure  3.  (A)  A  slice  of  the  potential  energy  surface  along  the  diagonal  of  the  normal  mode  coordinate,  Q,  as  each  substituted  pentacene  monomer  in  the 
dimer  relaxes  along  the  1435  cm-1  mode.  The  values  at  each  point  are  calculated  at  the  8o8e  level  of  theory  (see  the  Supporting  Information  for  details). 
The  ME  state  is  higher  than  the  LE  states  at  Q.=  (0,0),  and  the  avoided  crossing  is  Q_=  (0.33,0.33).  (B)  Percent  diabatic  state  composition  of  the  three 
lowest  excited  adiabatic  states  of  BP-37  as  a  function  of  (^demonstrating  strong  mixing  at  Q_=  (0.33,0.33)  and  the  change  in  character  of  adiabatic  states 
S[  and  S3  as  we  move  from  Q_=  (0,0)  to  Q_=  (1.3, 1.3).  (C)  Representation  of  the  1435  cm-1  mode  in  pentacene  and  the  postulated  structure  of  T,  in  a 
pentacene  monomer,  achieved  through  the  avoided  crossing. 


Table  3.  Energies  (in  eV)  and  Diabatic  Contributions  (in  %)  of  the  Eight  Lowest  Adiabatic  Singlet  States”  Calculated  at  the  4o4e 
Level  of  Theory 
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“See  Figure  2  for  representation  of  the  diabats. 


concomitant  with  the  elongation  of  the  other  bonds.  The  bond¬ 
stretching  and  contraction  suggests  stabilization  of  a  biradical 
structure  with  triplet  character  localized  in  the  central  pentacene 
ring  (animation  of  Q_  available  as  Supporting  Information).  This 
biradical  character  results  in  rapid  lowering  of  the  ME  state 
energy,  resulting  in  the  near-degeneracy  (avoided  crossing) 
necessary  for  strong  direct  mixing  of  the  LE  and  ME  states.54-58 
It  must  be  noted  that  the  structure  on  the  right  is  not  a  minimum 
energy  structure,  but  one  that  represents  the  bond  character 
along  the  important  mode  in  iSF.  It  is  clear  from  the  analysis 
above  that  the  rate  of  iSF  will  be  intimately  tied  to  the  direct 
coupling  matrix  elements.  This  coupling  can  be  tuned  by  the 
extent  of  orbital  overlap,  in  particularly  via  geometric  changes. 
For  instance,  as  the  monomers  rotate  around  the  central  bond 
and  approach  orthogonality,  we  expect  coupling  to  reach  a 
minimum.  We  consider  a  rotation  of  the  dihedral  angle  from  37° 
to  57°  and  find  a  2-fold  decrease  in  the  direct  coupling  matrix 
element  at  the  avoided  crossing  region  (~3  meV  vs  ~1.5  meV). 


Assuming  a  golden-rule  like  expression  for  the  rate  of  singlet 
fission  where  the  rate  is  proportional  to  the  square  of  the  direct 
coupling  strength,4  we  predict  a  decrease  in  the  iSF  rate  by  a 
factor  of  4;  a  key  result  that  will  be  compared  to  the  predicted 
CT-mediated  value  in  the  next  section,  and  confirmed  by 
experimental  results  below. 

■  QUANTIFYING  THE  ROLE  OF  CHARGE  TRANSFER 
STATES 

While  the  identification  of  an  avoided  crossing  is  essential  to 
make  the  direct  coupling  mechanism  viable,  this  does  not  in  any 
way  rule  out  a  contribution  from  the  CT-mediated  super¬ 
exchange  mechanism.  In  this  section,  we  investigate  the  role  of 
CT  states  and  we  quantify  the  change  in  iSF  rate  associated  with 
the  CT-mediated  mechanism  as  the  dihedral  angle  is  varied. 

As  previously  mentioned,  we  find  that  accurate  quasi-diabatic 
CT  states  cannot  be  obtained  at  the  8o8e  level  of  theory.  In  order 
to  quantify  the  extent  of  CT  involvement  we  move  to  a  reduced, 
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Figure  4.  Energies  of  the  4o4e  diabatic  states  as  a  function  of  (A)  intermonomer  C— C  bond  length  in  BP-37  and  (B)  the  dihedral  angle  shown  in  Figure 
IB.  (C)  CT-mediated  superexchange  VeS  as  defined  in  the  text  as  a  function  of  the  bipentacene  dihedral  angle. 


four-orbital/four-electron  (4o4e)  active  space  that  includes  only 
the  HOMO  and  LUMO  on  each  pentacene  for  BP-37.  Table  3 
lists  the  adiabatic  state  energies  and  the  corresponding  diabatic 
state  contributions  from  4o4e  calculations  and  illustrates  the 
minimal  mixing  observed  between  different  “types”  of  diabatic 
states. 

The  adiabatic  states  largely  correspond  to  a  single  type  of 
diabatic  state,  with  the  exception  of  S0  and  S4;  where  the  inherent 
biradical  character  of  BP-37  results  in  low  energy  doubly  excited 
states,  dg  and  gd,  mixing  with  the  gg  diabat.  In  the  less  accurate 
4o4e  level  of  theory,  Si  is  largely  composed  of  ft  character,  while 
the  other  states  exhibit  linear  combinations  of  similar  types  of 
diabats  (eg  +  ge,  dg  +  gd,  or  ac  +  ca).  This  is  in  agreement  with  the 
8o8e  calculations  where  we  see  no  CT  character  in  states  within 
~1  eV  of  the  low-lying  singlet  excited  states.  We  also  find  that 
calculations  at  the  4o4e  level  underestimate  energies,  particularly 
the  ME  state,  but  are  able  to  reproduce  the  qualitative  picture  of 
high-lying  CT  states  and  very  little  diabatic  mixing. 

Previous  work  investigating  the  mechanism  ofxSF  unfavorably 
compared  the  direct  coupling  mechanism  against  superexchange 
mediated  by  high-lying  CT  states.2' -  1  Here,  the  exper¬ 
imentally  observed  0.76  ps  iSF  time  constant  in  BP-37,1  the 
large  energy  separation  (~1  eV)  of  the  CT  states,  and  the 
relatively  small  energy  separation  between  the  ME  and  EE  states 
suggest  a  mechanism  that  may  be  independent  of  CT  mediation. 
Although  the  magnitude  of  the  coupling  between  the  ME/LE 
states  and  the  CT  states  in  BP-37  is  comparable  to  values 
reported  for  pentacene  and  tetracene  dimers  that  undergo 
xSF,  the  energetic  separation  of  the  C 1  states  tor 

iSF  in  BP-37  is  significantly  larger.12’31’32’34  This  difference  is 
expected  since  the  dimer  lacks  the  HOMO/LUMO  energy 
bands  found  in  crystalline  pentacene,  in  addition  to  the  longer 
intermonomer  distance  (monomer  center-of-mass  separated  by 
13.6  A  in  BP-37  vs  5  A  in  crystalline  pentacene)  that  leads  to  a 
reduction  in  stabilizing  Coulomb  binding  energy. 

(fflHlac)(adHleg)  +  (fflHlca)(calHleg) 

'  "  ”  Ect  -En  +  Ect  -  Eeg  (i) 

Unlike  the  direct  coupling  mechanism,  the  CT-mediated 
mechanism  is  second-order  in  nature,  and  depends  on  two 
factors:  the  coupling  matrix  elements  and  energy  differences 


between  states.  The  effective  indirect  coupling  is  defined  by  eq  1 
and  likewise  for  the  ge  diabat. 

Because  CT  mediation  is  an  inherently  nonlocal  phenomenon, 
it  is  expected  that  CT  energies  and  the  indirect  coupling  may 
depend  strongly  on  the  intermonomer  geometry  of  the 
covalently  linked  pentacene  units.  We  examine  two  degrees  of 
freedom  that  affect  the  interactions  between  the  pentacene 
chromophores:  the  length  of  the  covalent  bond  between  the  two 
monomers  and  the  dihedral  angle  6.  We  find  that  CT  state 
energies  are  at  a  minimum  near  the  reference  C— C  bond  length 
of  1.48  A,  and  increase  only  slightly  in  energy  as  the  a  bond 
length  decreases  (Figure  4A).  Similarly,  Figure  4B  shows  that  the 
energies  of  the  CT  states  remain  high  in  energy  and  fairly 
constant  relative  to  the  lower  lying  excited  states  as  a  function  of 
dihedral  angle. 

Although  CT  states  are  substantially  higher  in  energy,  this  does 
not  necessarily  preclude  an  efficient  CT-mediated  mechanism. ,0 
A  better  method  to  establish  the  role  of  CT  state  mediation  is  to 
examine  the  effective  coupling  as  a  function  of  chemical 
structure,  and  compare  the  predictions  of  this  mechanism  to 
the  experimentally  observed  iSF  rates.  Drawing  parallels  to  our 
analysis  of  the  direct  mechanism,  we  show  the  indirect  coupling 
as  a  function  of  dihedral  angle  in  Figure  4C.  Further,  the  indirect 
coupling  decreases  by  a  factor  of  ~6  as  the  dihedral  angle  changes 
from  37°  to  57°,  while  the  CT  state  energies  remain  roughly 
constant  relative  to  the  ME  and  LE  sates.  Assuming,  again,  a 
golden-rule  rate  expression,  this  would  predict  slower  iSF  by  a 
factor  of  ~36  in  contrast  to  the  factor  of  4  predicted  by  the  direct 
coupling  pathway.  We  note  that,  despite  the  high-level 
multireference  methods  used  here,  we  cannot  directly  compare 
the  magnitude  of  coupling  for  the  direct  and  CT-mediated 
mechanisms,  rather  we  rely  on  trends  in  coupling  to  predict 
changes  in  iSF  rates.  We  now  turn  to  experimental  evidence  to 
distinguish  between  these  two  different  scaling  relationships  and 
to  identify  the  dominant  mechanism. 

■  EXPERIMENTAL  RESULTS 

In  order  to  experimentally  investigate  how  the  excited  state 
dynamics  are  affected  by  modifying  interplanar  angle  between 
the  two  substituted  pentacene  monomers,  we  synthesize 
compounds  BP-HPh  ( 6  =  42°)  and  BP-2Ph  (where  0  is 
analogous  to  BP-57),  where  the  dihedral  angle  twist  is  imposed 
by  steric  hindrance  from  the  R  groups  shown  in  Figure  IB.  Since 
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Figure  5.  (A— C)  Raw  transient  absorption  (AA)  data  for  dilute  solutions  of  BP-2H,  BP-HPh,  and  BP-2Ph  in  chloroform  plotted  as  a  function  of 
wavelength  and  time.  Regions  of  large  pump  scatter  at  ~600  nm  have  been  removed  for  clarity.  (D— F)  Corresponding  singlet  and  triplet  spectral 
deconvolution  for  BP-2H,  BP-HPh,  and  BP-2Ph,  respectively,  solved  by  singular  value  decomposition  and  global  analysis. 


the  phenyl  groups  are  twisted  orthogonal  to  the  pentacene  core, 
they  do  not  have  significant  contributions  to  the  electronic 
structure,  as  evidenced  in  the  UV— vis  spectrum  where  their 
addition  does  not  shift  the  onset  of  absorption  (i.e.,  the  phenyl 
groups  do  not  participate  in  delocalization  with  pentacene).  The 
transient  absorption  (TA)  spectra  (Figure  5)  of  BP-HPh  and  BP- 
2Ph  show  qualitatively  similar  features  to  the  parent  BP-2H 
(which  is  analogous  to  BP-37),  and  the  iSF  rates  and  yields  can  be 
similarly  obtained.  The  spectrally  deconvoluted  singlet  and 
triplet  transient  features  are  evaluated  using  global  analysis 
methods,6"  and  all  three  exhibit  a  broad  singlet  excited  state 
transition  that  spectrally  overlaps  with  a  narrower  triplet 
absorption  feature  at  ~520  nm.  A  triplet  absorption  feature  on 
the  red  side  of  the  ground  state  bleach  is  prominent  in  BP-2H  but 
decreases  in  intensity  as  the  molecule  becomes  more  twisted.  We 
have  seen  a  similar  effect  in  bipentacene  molecules  with 
conjugated  spacers  and  believe  that  this  trend  can  be  attributed 
to  coupling  of  the  Tj  — >  T„  triplet  transition  dipole  moments, 
oriented  along  the  long  axis  of  the  molecule  when  the  pentacenes 
are  both  aligned  (planar)  and  in  close  proximity.6^ 

As  in  BP-2H,  we  quantify  the  iSF  yield  from  the  ground  state 
bleach  signal  in  the  deconvoluted  spectral  transients  for  the 
singlet  and  triplet.15’18'64  This  method  is  valid  since  the  entirety 
of  the  molecule  is  bleached  when  populated  by  one  singlet  or  two 
triplets.15’1  We  establish  that  increased  twisting  of  the  molecule 
does  not  violate  this  condition  by  measuring  each  sample  under 
the  identical  excitation  conditions  and  verifying  that  the  quantity 
AA/(e(l  —  Tj)  is  conserved  for  the  ground  state  bleach,  where  e 
is  the  extinction  coefficient  at  the  absorption  maximum  and  T  is 
the  transmission  at  the  pump  wavelength.  In  other  words,  the 
singlet  is  fully  delocalized  on  the  whole  molecule. 

We  observe  that  the  area  of  the  ground  state  bleach,  when 
corrected  for  overlap  with  induced  absorption  signals,  is 
unchanged  after  singlet  fission.  Thus,  all  initially  populated 
states  are  converted  to  triplet  pairs  and  a  quantitative  iSF  yield  of 
~200%  is  observed.  This  high  yield  is  consistent  with  a  simple 
rate  analysis,  since  singlet  fission  is  orders  of  magnitude  faster 
than  competing  deactivation  processes  such  as  the  ~12.3  ns 
fluorescence  decay  in  TIPS-pentacene  monomer.6  .  While 
qualitatively  similar  dynamics  are  observed  in  all  three 
compounds,  quantitative  differences  in  the  rates  of  singlet  fission 
and  triplet  recombination  are  observed.  In  Figure  6,  we  show  the 


Figure  6.  Single  wavelength  transient  absorption  kinetics  monitoring 
the  peak  in  the  triplet  excited  state  absorption  spectrum  in  chloroform, 
between  0.1  and  0.2  mW  and  less  than  100  fiM  for  all  samples.  While 
both  iSF  and  the  corresponding  triplet  decay  rates  slow  with  increased 
interplanar  angle,  there  is  a  larger  change  in  the  rate  of  triplet  decay. 


TAS  kinetics  associated  with  the  peak  of  the  triplet  absorption  in 
chloroform  at  517,  516  and  528  nm  for  BP-2H,  BP-HPh  and  BP- 
2PH,  respectively.  It  should  be  noted  that  there  is  some  overlap 
with  induced  absorption  from  the  singlet  state,  which  gives  an 
instantaneous  rise  at  time  zero.  Global  fitting  of  the  full  data  set  is 
used  to  achieve  high  accuracy  of  the  extracted  time  constants. 

The  ultrafast  time  constant  for  iSF  monotonically  increases  as 
the  dihedral  angle  increases,  from  0.76  ps  to  1.69  ps  to  3.38  ps. 
The  iSF  rate  constant  changes  from  BP-2H  (least  twisted)  to  BP- 
2Ph  (most  twisted)  by  a  factor  of  ~4.5,  very  close  to  the  factor  of 
4  predicted  by  the  direct  coupling  mechanism  comparing  BP-37 
to  BP-57,  and  an  order  of  magnitude  different  from  the  slow¬ 
down  in  rate  predicted  by  the  CT-mediated  superexchange 
mechanism.  We  note  that  the  rate  of  triplet  recombination  is 
similarly  affected,  with  the  time  constant  monotonically 
increasing  from  450  ps  to  1600  ps  to  5200  ps,  as  the  dimer 
becomes  more  twisted.  We  find  that  the  change  in  electronic 
coupling  strength  for  the  recombination  of  two  triplets  appears 
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to  be  more  dramatic  than  for  iSF,  with  recombination  rates 
scaling  even  more  strongly  than  iSF  rates  with  the  twist-induced 
reduction  in  coupling.  Since  fast  iSF  but  slow  triplet 
recombination  are  desired,  this  observation  suggests  a  viable 
way  to  optimize  the  overall  excited  state  lifetime,  an  important 
result  for  technological  applications. 

■  EFFECT  OF  SOLVENT  PROPERTIES  ON  ISF  RATES 

Finally,  using  both  theory  and  experiments,  we  establish  that 
solvent  properties  do  not  significantly  and  systematically  change 
the  rate  of  iSF,  suggesting  that  iSF  proceeds  primarily  through  a 
direct  coupling  mechanism  in  the  bipentacenes  studied  here.  We 
calculate  diabatic  excited  state  energies  in  solution  by  employing 
the  polarizable  continuum  model  (PCM)  to  simulate  solvent 
effects.4' '  Even  for  solvents  as  polar  as  water,  we  observe  only  a 
minimal  50  meV  shift  in  CT  state  energies  relative  to  all  other 
states  (further  details  provided  in  the  Supporting  Information). 
Experimental  studies  reveal  a  weak  and  complex  dependence  of 
the  electronic  energy  levels  on  the  properties  (in  particular 
orientation  polarizability)  of  the  solvent,  consistent  with  the  idea 
that  direct  coupling  mediates  iSF  and  CT  states  are  not 
necessarily  involved.  While  small  changes  in  the  energy  of  the  S0 
— >  S,,  Si  — ►  SN,  and  T,  — »  T,  transitions  are  observed,  there  is  a 
nonmonotonic  dependence  on  the  solvent  dielectric  constant, 
index  of  refraction,  or  orientation  polarizability  (Supporting 
Information).  The  singlet  fission  rate  constants  change  by  less 
than  a  factor  of  2  over  the  entire  range  of  solvents  used  (including 
hexane,  xylenes,  chloroform,  THF,  acetone),  which  is  consistent 
with  a  dominant  mechanism  that  does  not  involve  CT  states  to 
mediate  iSF.  We  note  that  other  pentacene  dimers  with  a 
markedly  different  connectivity  were  recently  shown  to  exhibit 
weakly  solvent-dependent  iSF  dynamics.1  ,lx  That  result 
suggests  that  molecules  with  varying  degrees  of  conjugation 
and  through-space  interactions  may  adopt  intramolecular 
conformations  capable  of  accessing  lower  lying  CT  states. 

In  summary,  we  find  a  novel  mechanism  for  ultrafast  iSF  in 
pentacene  dimers:  direct  coupling  of  the  singlet  to  multiexciton 
state  facilitated  by  vibronic  coupling  induced  degeneracy. 
Furthermore,  we  rule  out  dominant  contributions  of  a  CT- 
mediated  mechanism,  unlike  in  the  case  of  xSF.  Experimental 
characterization  of  iSF  rates  of  molecules  with  different  dihedral 
angles  confirms  that  the  direct  coupling  mechanism  accurately 
reproduces  measured  trends  and  further  raises  the  interesting 
point  that  slowing  the  rate  of  iSF  in  highly  twisted  dimers  is 
accompanied  by  dramatically  slower  rates  of  triplet  pair 
recombination.  Establishing  the  mechanism  for  iSF  in  these 
materials  is  the  first  step  toward  designing  new  materials  with 
optimized  excited  state  properties  for  high  efficiency  organic 
solar  cells. 

■  EXPERIMENTAL  SECTION 

Calculations.  All  geometry  optimizations  described  in  the 
manuscript  are  performed  using  DFT-B3LYP  with  a  6-31G* 
basis  set.  Following  earlier  work  on  xSF  in  pentacene," "  excited 
state  energies  for  the  dimer  are  obtained  from  multiconfiguration 
quasi-degenerate  perturbation  theory  (XMCQDPT)  calcula¬ 
tions1"  with  a  0.02  au  intruder  state  avoidance  shift4  ’  using  the 
SBKJC  pseudopotential  with  its  corresponding  double-f  basis 
set  1  plus  a  d-polarization  function  taken  from  the  standard  6- 
31G*  basis  set  for  carbon.1'  We  find  that  an  8o8e  active  space  is 
necessary  to  yield  accurate  excited  state  energies  for  the  dimer,  a 
conclusion  similarly  reached  for  related  acene  dimers.1  ,4’'44 
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Diabatization  of  the  eight  (four)  lowest  lying  states  is  performed 
at  the  4o4e  (8o8e)  level  employing  the  4-fold  approach  of 
Nakamura  and  Truhlar  to  obtain  our  diabatic  states.68’69 
Reference  orbitals  are  generated  from  a  Hartree— Fock 
calculation  at  large  intermonomer  distances  which  allows  for 
each  orbital  to  be  unambiguously  associated  with  only  one 
monomer.  To  investigate  the  effects  of  solvent  polarity, 
diabatization  calculations  utilizing  the  polarizable  continuum 
model  (PCM)  are  performed  at  the  CASSCF  level  using  a 
minimal  active  space  so  as  to  explicitly  construct  charge  transfer 
diabats.48'65  A  state  specific  shift,  defined  as  E(XMCQDPT)  — 
E(CASSCF)  in  vacuo,  is  applied  to  the  resulting  energies  to 
account  for  the  absent  dynamical  correlation.  All  calculations  are 
performed  using  the  electronic  structure  package  GAMESS- 
US.70’71 

Transient  Absorption  Spectroscopy.  Transient  absorp¬ 
tion  spectroscopy  was  performed  using  a  commercial  Ti:sapphire 
laser  system  (SpectraPhysics)  operating  at  a  repetition  rate  of  1 
kHz.  A  commercial  optical  parametric  amplifier  (LightConver- 
sion)  was  used  to  generate  resonant  pump  pulses  with  ~100  fs 
time  duration.  Supercontinuum  probe  light  was  generated  by 
focusing  a  small  portion  of  the  800  nm  fundamental  into  a 
sapphire  disk.  The  probe  light  was  split  into  signal  and  reference 
beams,  which  were  detected  on  a  shot-by-shot  basis  by  a  fiber- 
coupled  silicon  (visible)  or  InGaAs  (infrared)  diode  array.  The 
pump— probe  delay  was  controlled  by  a  mechanical  delay  stage. 
Excitation  fluence  in  each  measurement  was  approximately  50 
^J/cm2.  Method  details  and  results  for  additional  ultrafast 
experiments,  including  triplet  photosensitization  and  photo¬ 
luminescence  upconversion  measurements,  can  be  found  in  the 
Supporting  Information. 
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■  ABBREVIATIONS 

SF,  singlet  fission;  xSF,  intermolecular  singlet  fission;  CT,  charge 
transfer;  iSF,  intramolecular  singlet  fission;  LE,  local  excitonic; 
ME,  multiexcitonic;  CASSCF,  complete  active  space  self- 
consistent  field  theory;  XMCQDPT,  extended  multiconfigura¬ 
tion  quasi-degenerate  perturbation  theory;  PES,  potential  energy 
surface;  ORMAS-SF,  occupation  restricted  multiple  active  space 
spin  flip;  8o8e,  eight  orbitals/8  electrons;  4o4e,  four  orbitals/four 
electrons;  PCM,  polarizable  continuum  model;  THF,  tetrahy- 
drofuran 
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